Paradoxical non-linear response of a Brownian particle 
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We consider a Brownian particle in a "meandering" periodic potential when the ambient temper- 
ature is a periodically or stochastically varying function of time. Though far from equilibrium, the 
linear response of the particle to an external static force is exactly the same as in the equilibrium 
case, i.e. for constant temperature. Even more surprising is the non-Unear response: the particle 
slows down and then even starts to move in the direction opposite to the applied force. 
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The second law of thermodynamics requires that the 
response of a system at thermal equilibrium to a static 
external force is a motion in the direction of that force. 
In this Letter we address the geometrically constrained 
motion of a Brownian particle with the remarkable prop- 
erty that, though being far from equihbrium, its linear 
response to a static force F is exactly the same as in the 
equihbrium case. In particular, the system is at rest in 
the absence of the external perturbation F and for suffi- 
ciently small (positive or negative) F it moves, as naively 
expected, in the same direction as the static force. An 
even more remarkable, genuine non-equilibrium behavior 
arises in the non-linear response regime: upon increas- 
ing F to moderately large (positive or negative) values, 
the particle slows down and then even starts to move in 
the direction opposite to the applied force F ! Note that 
such a response behavior is fundamentally different from 
so-called absolute negative conductance or mobility [l|, 
characterized by linear and non-linear response proper- 
ties that are just reciprocal to those described above, 
and from differential negative mobility [3], where the 
transport speed decreases with increasing force but never 
points in the direction opposite to that force. Accord- 
ingly, also the underlying physics and, in particular, the 
considered systems are entirely different. Furthermore, 
our present case has nothing to do with a so-called ratchet 
effect Q, where, by definition, the velocity of the Brow- 
nian particle is finite even in the absence of an external 
force, thus keeping the same sign in an entire interval of 
positive and negative forces F^ and moreover increasing 
monotonically as a function of F . 

The above announced response of a Brownian particle 
to an external force requires the interplay of two indis- 
pensable ingredients: a non- linear dynamics and a source 
of disequilibrium. Both will be chosen in a way that is 
conceptually very simple and experimentally easy to re- 
alize. As for the non-linear dynamics, we consider an 
overdamped Brownian particle, whose motion is confined 
to a quasi-one-dimensional meandering path in the x-y- 
plane, and which is subjected to an externally applied 
homogeneous force along the a:-axis, see Fig. 1. Experi- 
mentally, such a meandering state space can be realized 
e.g. by means of light forces JH or morphologically via 
lithographic etching methods |^ , while the external force 
can be generated e.g. by electrical or magnetic fields. 



or by gravitation. The second main ingredient of our 
model is a source of disequilibrium. Since the heydays of 
Carnot's machines, one of the simplest examples in this 
respect is a temperature that periodically switches be- 
tween two different values, see Fig. 2a. Apart from this 
"minimal" example we will also admit somewhat more 
general time dependent temperatures T(i), consisting of 
a "background" or "floor" temperature Tq , superimposed 
by short, but violent temperature "pulses", "spikes", or 
"outbursts" , occurring at random times ti and extending 
over characteristic time intervals r^, see Fig. 2b. In anal- 
ogy to the periodic case, the average interspike distance 
is denoted by 

t:={U+i-U), (1) 

and the relative duration of the temperature spikes by 



^ := {n)/T 



(2) 



Experimentally, such temperature variations may be re- 
alized either directly via heating e.g. with laser pulses 
or indirectly e.g. by pressure pulses. We also note that 
such outbursts may not necessarily be caused by ther- 
mal noise but may as well be considered as a simplified 
model for sudden mechanical "strikes" , "quakes" , "im- 
pacts", "eruptions" etc. With a neghgibly small "back- 
ground" temperature Tq, such a model may then describe 
a variety of experiments with even quite large "particles" 
subject to some kind of sporadic violent shaking. 

The quantity of central interest is the average particle 
current along the a;-axis 



lim 



x{t) 



t 



(3) 



For ergodicity reasons, no ensemble average on the right 
hand side is needed and due to the long time limit, initial 
conditions do not matter. For the symmetric example in 
Fig. lb it is clear that the current (i) is an odd function 
of the applied force F and in particular vanishes in the 
absence of the force. For more general cases like in Fig. 
la, the former property is lost but, as we will show later, 
one still has (i) = for F = 0. In other words, our 
system docs not exhibit a ratchet effect. 

Parametrizing by q the position of the particle along its 
meandering route (path length), its overdamped Brow- 
nian motion can be modeled by the Langevin-equation 



iH 



77 m = -r{q{t)) + v^k^m m , 



(4) 



where 77 denotes the damping coefhcient, fee Bohzmann's 
constant, and ^(i) a 5-correlated, unbiased Gaussian 
noise. The force —V'{q) in Q) represents the projec- 
tion of CxF in Fig. 1 along the meandering path at 
the position parametrized by q. Thus V{q) is propor- 
tional to F. Denoting by Lq the natural spatial period 
along the meandering path one can see from Fig. 1 that 
increasing the path length q by one period Lg is tan- 
tamount to increasing the x-component of the particle 
position by one period L^ and leaving the y-component 
unchanged. As a first consequence, we can infer that 
V{q + Lq) — V{q) - FLj. and hence 



V{q) 
G{q + Lq) 



F[G{q) 
G{q). 



qLx/Lq 



(5) 
(6) 



As a second consequence, it follows that the average ve- 
locity of actual interest © is related to the average ve- 
locity (q) := \iTat^oo{q{t) /t} of the auxiliary dynamics 
Q according to 



(i) = (q) L,x/Lg 



(7) 



Note that the periodic function G{q) from (jsj, lO is 
-F- independent. In Fig. 3, we exemplify its explicit con- 
struction for the piecewise linear meandering path from 
Fig. lb. Whence one can infer that 



G{q) 



-q{l-Lx/Lq) 
q{cosip + Lx/Lq) 



for -I2 
for 



< q < 

< q < h 



(8) 



G{q+ h +I2) = G{q) 
Lx ~ 2 {I2 ~ h cosif) 
Lq = 2{h + h) . 



(9) 
(10) 

(11) 



In order to study the linear response behavior of our 
system, the deterministic force —V'{q) in Q, being Lg- 
periodic according to ((SJ, 0, is expanded into a Fourier- 
series. Since each summand of this Fourier expansion is 
proportional to F (cf. Eq. lO) it follows that in lead- 
ing order F (linear response) the effect in (0} of each of 
those summands can be considered separately and their 
individual contributions to the current (q) can simply 
be added up. For symmetry reasons one can infer y 
that any summand which exhibits harmonic oscillations 
as a function of q does not lead to a systematic directed 
transport, i.e. does not contribute to (q). It remains 
the non-oscillating part of —V'{q)^ which according to 
Q is given by FL^/Lg, and hence gives rise to a cur- 
rent ri{q) = FLx/Lq after averaging out the noise in Q). 
Exploiting iQ we thus can infer that 



(i) 



F 



0{F^) 



(12) 



The same result can also be derived in a mathematically 
more rigorous way by means of a perturbation expansion 
in F of the Fokker-Planck equation equivalent to Q. 
Moreover, it can be extended beyond the so far consid- 
ered case that T{t) > for all t. The most remarkable 
property of the linear response behavior (|12|l is its inde- 
pendence of the temperature T(i): For strong tempera- 
ture variations, i.e. far from equilibrium, it is exactly the 
same as in the equilibrium case when the temperature is 
constant! 

Next we turn to the non-linear response for not too 
small forces F. The simplest case arises when F is so 
large or the "background temperature" Tq is so small 
that in the auxiliary dynamics Q-JHll the effect of the 
fluctuations ^{t) is negligible between successive temper- 
ature spikes, cf. Fig. 2. Furthermore, any temperature 
outburst is assumed to be so violent, that the distribution 
of an ensemble of Brownian particles will be completely 
randomized, i.e. practically uniformly distributed within 
any interval comparable to the spatial period Lg. More- 
over, during such an outburst, the particles will practi- 
cally not feel the details of the periodic part G{q) of the 
potential in (O but only its systematic "tilt" —qFL^/Lq. 
Hence, the average displacement during an outburst of 
duration r^ is 

Aqsin) = nFLj{Lq7j) . (13) 

Between two outbursts, the particles start with a uniform 
initial probability distribution and then evolve determin- 
istically according to Q) with ^(t) = 0. Focusing on the 
piecewise linear example from Fig. 3b and Eqs. l|S|l- Hll|) 
and on positive forces F, one readily finds that the aver- 
age net displacement during the time Ai^ := t^+i —ti — n 
between two temperature outbursts amounts to 

AqoiAt,) = [q - {h - v2At,fQ{l2 - V2At,)]/Lg 

- [ll - {h - viAt,fQ{h - viAti)]/Lg (14) 

where V2 '■= \F\/r] and vi := \F\ cosip/r] are the two de- 
terministic particle speeds associated with the two dif- 
ferent slopes of the piecewise linear potential V{q) in 
Fig. 3b and the Heaviside step functions Q account for 
the fact that the particles stop to move once they have 
reached a local potential minimum of V{q). The average 
displacement per temperature outburst then follows as 
{Aqs{Ti) + A(7o(Aii)). Dividing by the average interspike 
distance from IJT} yields {q). With ©,0, CSl) and tak- 
ing into account that (i) is obviously an odd function of 
F, we finally obtain 

''^^ , (Ago(tj+i -tj -Tj)) ' 
L„ \F\t 



(i) 






V 



(15) 



For random temperature bursts, the average (A(7o(ii+i — 
ti — Tj)) has to be evaluated according to their specific 
statistical properties. In many cases, one expects that the 
result will be well approximated by Aqo{{ti^i —tt—Ti)) = 
Ago([l-^]T),i.e. 



{x)^F- 



L, 



Aqo{[l-i)]T) 
\F\t 



(16) 



For periodic bursts, Eqs. H15() and H16|l are even strictly 
equivalent. 

For small F, one readily recovers the correct linear 
response behavior 1)12(1 from (|15|) . though this small- i^ 
regime was originally not included in our considerations 
above. Next we turn to the case that |F| is sufficiently 
large that the 0-functions in H15|l vanish. Focusing on 
the case li > I2 in Fig. 3, this amounts to the condition 
(1 - d)T\F\ cos If /rj > li, and (23 with (^, (^ then 
takes the form 



±Jrp 






2\F\t 



(17) 



For sufficiently short relative pulse durations •& we thus 
expect the existence of an interval of moderately large 
\F\-values with the property that the current (x) will be 
positive for negative F and vice versa! 

In Figs. 4 and 5, results of numerical simulations are 
presented, exhibiting a fairly good agreement with the 
approximation H15|l . In particular, both the predicted 
linear response and the paradoxical non-linear response 
behavior are confirmed. While for the filled circles the 
assumptions made in the derivation of the approximation 
(|15|l are apparently quite well satisfied, for the open cir- 
cles there are notable deviations. As mentioned above, 
for the periodic T{t) in Fig. 4 the approximations H15|) 
and 116() coincide, whereas for the stochastic T{t) in Fig. 
5 the more complicated expression (|15|l is clearly supe- 
rior. 

Finally, we turn to a heuristic discussion of the general 
conditions under which the above paradoxical non-linear 
response behavior is expected. Taking into account the 
obvious constraints < if < 7r/2, Z2/2 — li cost/? > in 
Fig. 3 and Zi > I2, {I — d)T\F\ cos (f/r] > li [see above Eq. 
(|17|l ] one can infer that a necessary condition to get a cur- 
rent opposite to F in H17|) is ?? < 1/3. On the other hand, 
the shorter the relative duration i? of the temperature 
bursts in (|17|l are, the more pronounced the effect will 
be. In particular, for i? ^ [i.e. 5-spikes in Fig. 2] the 
current (i) in (|17|l remains opposite to F for arbitrarily 
large |_F|. So far we have assumed, that each temperature 
burst results in an almost uniform particle randomization 
within any period Lq. In the opposite case, i.e. for not 
so violent bursts, the distribution will remain notably 



peaked around the local potential minima of V{q), lead- 
ing to a reduced displacement opposite to F during the 
subsequent relaxation period till the next burst. In other 
words, the effect of interest diminishes with decreasing 
intensity of the temperature bursts, and obviously disap- 
pears altogether when the temperature variations in Fig. 
2 become negligible. Similarly, one sees that our assump- 
tion of a negligibly small "background" temperature Tq, 
i.e. a deterministic relaxation dynamics in between suc- 
cessive bursts, is optimal in the sense that with increasing 
the noise intensity Tq, the effect of interest decreases and 
finally disappears. These predictions are confirmed by 
the open circles in Fig. 4, whose deviations from the 
theoretical solid line can be traced back to the increased 
To-value for small F and to the reduced T^igh-'^shie for 
large F. Finally, going over from the piecewise linear ex- 
ample in Fig. lb to the more general case in Fig. la leads 
to an increased deterministic relaxation time towards the 
local potential minima of the associated auxiliary poten- 
tials V{q) and hence again to a reduction of the effect 
of interest. Notwithstanding, it is clear that a current 
opposite to an applied force of suitable magnitude will 
arise whenever the outbursts are sufficiently short and 
violent, Tq is sufficiently small, and if the meandering 
path is such that a motion along this path with constant 
velocity will result in a back-and-forth motion when pro- 
jected along the x-axis, whose "backward-segments" are 
of longer duration than the "forward-segments" . Sim- 
ilarly, when going over from the quasi-one-dimensional 
paths in Fig. 1 (zero "width") to the corresponding me- 
andering structures with a finite "width" , we numerically 
verified (not shown) that the basic qualitative response 
behavior in Fig. 4 remains unchanged upon extending 
the width of the black lines in Fig. 3a at least as long 
as the corners of the structure indicated by the encircled 
numbers 1 and 4 in Fig. 3a did not yet merge (corre- 
sponding to a maximal width of about 0.182 /2). 
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FIG. 1: Examples of quasi-one-dimensional meandering 
paths in the x-y-p\a,ne (e.g. generated by a potential which 
is zero in the black and infinity in the white domains). They 
are required to be periodic (with period Lx) along the x-axis 
and confined along the y-axis. A further essential feature is 
the existence of sharp bends such that a steady motion along 
the meandering path results in a back-and-forth motion when 
projected along the a::-axis. In general, no further symmetry 
is required, as exemplified with (a). The piecewise linear ex- 
ample (b) of high symmetry serves to simplify the discussion 
of the basic physical mechanisms. A homogeneous, static bias 
force F is applied along the a;-direction. 
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FIG. 2: Examples of time-dependent temperatures T{t). (a) 
T{t) is periodically switching between Thigh and Tq. The total 
period is r and the duration of the high- and low-temperature 
segments j?r and (1 — ■!?)r, respectively, with i9 G (0, 1). (b) 
T{t) consist of a "background" To and random "spikes" at 
times ti of duration Ti. Their average distance and duration 
are characterized according to Eqs. O, lO- 
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FIG. 3; (a) Same example as in Fig. lb, parametrized by the 
angle <fi and the lengths h, I2. (b) The corresponding one- 
dimensional potential V{q) as described above Eq. Q, i.e., 
V'iq) = ~F for q € [-^2,0] or g £ [h,h + h] and V'{q) = 
Fcosip for q € [—h — hj—h] or q € [0,/i]. The encircled 
numbers mark corresponding positions in (a) and (b). 
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FIG. 4: Average particle current l^l versus static force F 
for the meandering path in Fig. 3a with h — 3, h ~ 2, 
if = 75°, and a time-periodic temperature according to Fig. 
2a with r = 10, and i9 = 0.1. Filled circles: To = 0.001, 
Thigh ~ 50. Open circles; To = 0.01, Thigh = 5. Shown are 
results from numerical simulations of CT- lllll in dimensionless 
units with 7] = ks ~ 1- The numerical uncertainty is about 
the symbol size. Solid line: The theoretical approximation 
1151 1. coinciding with II6II since T{t) is periodic. Dotted line: 
guide for the eye. 
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FIG. 5: Filled circles: Same as in Fig. 4 with the exception 
that the spiking times ti are not periodic [as in Fig. 2a] but 
rather randomly sampled [as in Fig. 2b] according to a Pois- 
son process with average interspike distance t — 10 [cf. 10] . 
The relative duration Ti = 0.07 = i9 is the same for all spikes i 
[cf. ||5J]. Solid line: Theoretical approximation 11511 . Dashed 
line: Simplified approximation H16|l . 



